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基于 MATLAB 的 多 元 非 线性 回归 模型 


董 大 校 
( 临 沧 师范 高 等 专科 学 校 ,云南 临 沧 6770000) 


摘 要 : MATLAB 是 源 于 和 矩阵 运算 的 一 种 高 度 集 成 的 计算 机 语言 。 它 提供 了 强大 的 科学 运算 .灵活 的 
程序 设计 流程 高 质量 的 图 形 可 视 化 与 界面 设计 ,便捷 的 与 其 他 程序 和 语言 接口 的 功能 。 文 章 充分 利用 
MATLAB 统计 工具 箱 的 优势 ,通过 程序 的 实现 ,对 多 元 非 线 性 回归 模型 的 未 知 参数 的 估计 方法 以 及 对 佑 
计 后 的 模型 预报 做 出 研究 ,并 以 实例 验证 了 该 方法 的 有 效 性 。 
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1 预备 知识 


非 线性 回归 最 小 二 乘法 拟 合 的 基本 原理 '"。 
对 给 定数 据 (x;,y;) (i = 0,1,…,m) ,在 取 定 的 函数 类 中 , 求 p(x) e ,使 误差 7 = p(x) 一 ii 
= 0,1,…,m) 的 平方 和 最 小 , 即 


= > [p(x;) -yi 
最 小 ,从 几何 意义 上 讲 , 就 是 寻求 与 给 定点 (x yi = 0,1,…,m) 的 距离 平方 和 为 最 小 的 曲线 y = 
p(x)( 图 1) ,函数 p(x) 称 为 拟 合 困 数 或 最 小 二 乘 解 , 求 拟 合 男 数 P(x) 的 方法 称 为 曲线 拟 合 的 最 小 二 乘 
法 。 


2 MATLAB 非 线 性 曲线 拟 合 命令 介绍 


2. 1 nlinfit 函数 王 

用 nlinfit 函数 进行 非 线性 最 小 二 乘 数据 拟 合 。 该 天 数 使 用 高 斯 一 牛顿 算法 ,调用 格式 如 下 : 

全 beta =njinfit(X,y,fun,beta0) 用 最 小 二 乘法 估计 非 线 性 函数 系数 。y 为 响应 值 ( 因 变 量 ) 矢量 。 
一 般 地 ,为 自 变量 值 组 成 的 设计 和 矩阵 ,每 一 行 对 应 与 y 中 的 一 发 个 值 ,但 是 ,X 可 以 是 fun 参数 能 接受 的 
任何 数组 。fun 参数 为 一 函数 ,该 函数 具有 下 面 的 形式 

yhat = myfun( beta ,X ) 

其 中 beta 为 系数 矢量 ,X 为 设计 矩阵。fun 为 参数 返回 一 个 拟 合 y 值 的 yhat 矢量 。beta0 为 包含 系 

数 初始 值 的 矢量 。 
全 |[ beta,r,J] =nlinfit(X,y,fun,beta0) 返回 拟 合 系数 (beta) 、 残 差 (r) .和 和 雅 可 比 矩 阵 J， 这 些 参数 
可 以 用 于 nlintool 函数 ,生成 预 页 测 值 的 误差 仿 计 ; 或 用 于 nlparci 枉 数 生成 系数 的 误差 估计 。 
2. 2 nlintool 函数 
使 用 nlintool 函数 可 以 对 数据 进行 非 线性 方程 拟 合并 交互 图 形 显示 ,其 调用 格式 如 下 : 
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@nlinfit(x,y,fun,beta0) 为 一 预测 图 , 它 提供 数据 (x,y) 的 非 线 性 曲线 拟 合 。 它 用 两 条 红色 曲线 
来 表示 预测 值 的 95% 置信 区 间 。beta0 为 一 矢量 ,包含 参数 的 初 值 。fun 参数 为 一 函数 ,该 函数 具有 yhat 
= myfun( beta ,X) 的 形式 : 

其 中 beta 为 系数 矢量 ,X 为 设计 和 矩阵 。 

nlintool 函数 显示 一 个 图 形 “ 和 天 量 ” ,其 中 的 每 个 图 形 对 应 于 输入 矩阵 x 的 每 一 列 。 响 应 变量 y 为 一 
列 矢量 ,与 x 的 行 数 相 匹 配 。 

@nlintool(x,y,FUN ,beta0 ,aljpha) ”用 图 形 表示 预测 值 的 100(1 -alpha)% 置信 区 间 。 当 x 为 矩阵 
时 ,将 为 每 一 列 单独 生成 图 形 ,应 变量 y 为 一 列 变量 ,对 应 于 x 的 行 。alpha( ) 的 默认 值 为 0.05 生成 
95% 置信 区 间 。 
2. 3 nlpredci 函数 

可 利用 njpredci 函数 计算 非 线性 模型 预测 值 的 置信 区 间 ,其 调用 格式 如 下 : 

@ypred = nlpredci( FUN ,inputs ,beta,r,J) 给 定 拟 合 参数 (beta) \ 残 差 (r) 和 和 雅 可 比 矩 阵 (J) , 返 
回 预测 响应 。 输 入 是 非 线 性 函数 中 独立 变量 的 数值 矩阵 。 
2.4 nlparci 函数 

可 利用 nlparci 函数 计算 非 线性 模型 中 参数 估计 值 的 置信 区 间 。 

@nlparci(beta,r,J) 给 定 拟 合 参数 (beta) 残 差 (r) .和 解 处 的 雅 可 比 矩 阵 (J) ,返回 非 线性 最 小 二 
乘 参数 估计 (beta) 的 95% 置信 区 间 ci。nlparci 函数 使 用 nlinfit 函数 的 输出 作为 输入 。 
2.5 regress 函数 

用 regress 盟 数 进行 多 元 线性 回归 。 

@|b,bint,r, rint, stats ] = regress(y,x,alpha) 给 出 bint 和 rint 的 100(1 -alpha)% 置信 区 间 。 
2. 6 stepwise 函数 

用 stepwise 旺 数 进行 逐步 回归 , 它 使 用 交互 环境 进行 分 析 。 其 调用 格式 如 下 : 

@stepwise(x,y,inmodel ,penter,premove) ”指定 模型 的 初始 状态 和 要 使 用 的 置信 区 间 。inmodel 为 
长 度 为 x 中 列 数 的 逻辑 矢量 。inmodel 指定 包含 在 初始 模型 中 的 自 变 量 。 默 认 时 不 包含 x 中 的 列 。pen- 
ter 指定 引入 目 变 量 的 最 大 Pp 值 ,默认 时 为 0. 05。premorve 指定 剔除 自 变量 的 最 小 p 值 ,默认 时 为 0. 10。 


3 实例 验证 


3.1 反应 动力 学 模型 

反应 动力 学 中 的 Hougen - Watson 模型 是 非 线性 模型 的 一 个 典型 例子 。 其 模型 如 下 : 

B! * *» — %3/B; 

l +pB, :x +pB;°*%, +B,*% 
其 中 ; B ,B,,… ,Bs 为 未 知 参 数 ,x ,x, 和 x 为 三 个 输入 变量 。 三 个 输入 为 氨 (hydrogen) ,n -成 烷 (n - 
pentane) 、 异 成 烷 (isopentane) 。 上 述 模型 显然 为 多 元 非 线 性 模型 。 文 件 reaction. mat 包含 了 反应 的 仿真 
数据 。 

通过 计算 可 知 ,反应 动力 学 Hougen - Watson 模型 中 的 5 个 未 知 参数 为 

B! = 1.2526,8, = 0. 0628,B, = 0.04 ,8, = 0. 1124 ,8;， = 1. 1914 


1.2526 . x, — x3/1. 1914 
最 终 模型 为 : rate = 1 + 0.0628 . x, +0.04 .x, +0.1124 .x, 
3. 2 财政 收入 预测 


财政 收入 与 国民 收入 .工农 业 总 产值 人 口 .就业 人 口 .固定 资产 投资 等 因素 有 关 。 采 用 逐步 回归 分 
析 构 造 预测 模型 。 以 财政 收入 作为 因 变量 y, 自 变量 为 国民 收入 (2 ) 、 工 业 总 产值 (x,) .农业 总 产值 


rate = 
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(%3) \ 总 人 口 (x) 就 业 人 口 (%s) 固定 资产 投资 (xs ) 。 其 样本 数 n = 30, 自 变量 数 p = 6 。 

首先 编写 非 线 性 匈 数 的 M 取 数 文件 model. m 代码 。 然 后 在 命令 窗口 输入 由 xi wx, xs3 xi xs 样本 数 
据 构 成 的 矩阵 x 和 财政 收入 样本 数据 y , 取 beta0 = [0.50 -0.03 -0.60 0.01 -0.02 0.35] ;用 交互 式 
非 线性 拟 合 图 形 工 具 来 显示 ,在 MATLAB 命令 窗口 中 输入 如 下 命令 : 

nlintool( X ,Y ,iodel beta0 ,0. 01 ) ;得 到 nlintool 分 析 窗 口 。 

由 此 得 到 财政 收入 与 各 因素 之 间 关 系 的 回归 模型 1 为 : 


y = 0.5245x, ~ 0. 0293x，- 0. 631x, + 0.0113x, ~ 0. 0231x，+ 0. 3652x。 (1) 
对 该 模型 进行 多 元 线性 回归 检验 ,输入 命令 : > > [b,bint,r,rint,stats] = regress(y,x) 
stats =0. 9805 301.0305 0 951.5530 


结果 分 析 :bint 为 各 系数 的 置信 区 间 。stats 矢量 是 值 分别 为 相关 系数 的 平方 .下 值 和 显著 性 概率 p。 
相关 系数 平方 值 R = 0. 9805 ,说 明 模 型 拟 合 程度 较 高 ,显著 性 概率 p = 0 ,小 于 0. 05 , 故 拒 绝 零 假 设 , 认 
为 模型 1 中 至 少 有 一 个 自 变 量 的 系数 不 为 零 ,因而 从 总 体 上 模型 1 是 有 意义 日 的 。 

结果 b 给 出 了 回归 系数 的 估计 值 ,检查 它们 的 置信 区 间 发 现 ,x, 系数 的 置信 区 间 包 含 了 零点 ,表明 
回归 变量 x, 对 变量 Y 的 影响 不 是 太 显著 ,因此 从 模型 中 移出 次 变量 ,对 该 模型 进行 逐步 回归 检验 ,在 
MATLAB 命令 窗口 输入 并 运行 命令 : > > stepwise (x,y) . 

除了 xi ,xs ,x6 外 ,其 他 县 变量 的 回归 系数 置信 区 间 都 包含 零点 ,保留 xi ,x, ,x6 ,将 其 他 变量 移出 回 
归 模 型 ,R 一 square =0. 979645 ,是 指 因 变量 的 97. 96% 可 由 模型 确定 , F = 417. 104 远 远 超过 了 下 检验 的 
临界 值 , = 0 远 小 于 显著 性 水 平 , 这 些 说 明 用 包含 % ,xs ,xs 的 线性 模型 来 描述 y 是 合适 的 oxi 的 回归 系 
数 为 0.382869 ,x 的 回归 系数 为 - 0. 496975 ,xs 的 回归 系数 为 0.408624 ,常数 项 由 Intercep 给 出 , 即 b = 
177. 07 ,利用 逐步 回归 得 到 的 回归 模型 2 为 : 

y = 0. 382869x; - 0. 49687Sx; + 0. 408624xs + 177. 07 (2) 

由 于 xi ,xs ,x6 对 y 具 有 较 明显 的 显著 性 ,用 更 为 精确 的 拟 合 曲线 来 拟 合 表示 x ,x3 ,xs 和 y 的 关系 。 
做 X] ,NX3 ,Ne 分 别 对 Y 的 散 点 图 ,具体 步骤 在 MATLAB 命令 窗口 输入 命令 : 

>>xl=x:,1);x3=x(:;3);x6o=x(:,6); > > subplot(1,3,1); plot(xl,y,’+’);> > subplot 

(1,3,2); plot(x3,y,’+);> > subplot(1,3,3); plot(x6,y,’+) 
再 做 经 变换 后 的 散 点 图 ,对 比 其 线性 效果 较 好 的 ,得 出 如 下 拟 合 曲线 : 
y =~71.391 +0.473x,;y = 1620. 46 - 6. 509x; + 0.01x3;y = 141.921 + 1. 158x。; 

在 MATLAB 命令 窗口 输入 如 下 命令 : 

> > subplot(1,3,1); > > plot(xl, —71.391 +0.473 *xl,’+); 

> > subplot(1,3,2); > > plot(x3 ,1620. 46 -6. 509 * x3 +0.01 * x3. 2 ,+ 

> > subplot(1,3,3); > > plot(x6 ,141.921 +1.1$8 * x6, +17 

对 应 的 尺 - Square 分 别 为 0. 951 ,0. 887 ,0. 944。 从 相关 系数 的 值 可 以 看 出 ,对 应 的 拟 合 方程 的 拟 合 
程度 较 高 。 所 以 ,用 拟 合 后 的 曲线 来 代替 原来 的 x, ,xs ,xs ,得 到 如 下 回归 模型 3: 

y = 0. 1811x, + 3. 2342x; ~ 0. 005x; + 0.4732x, — 597.4372 (3) 
作为 对 比 ,在 引入 交互 项 ,在 MATLAB 命令 窗口 输入 下 面 命令 : 

> >xl =x(:,1);x2 =x(:,2);x3=x(:,3);x7 =xl. *x2; x8 =xl. *x3; > > x=[xx7 x8];>> 
stepwise( x,y) 

加 入 交互 项 后 ,利用 逐步 回归 分 析 , 从 运算 结果 可 得 R - Square = 0. 98665 , = 461. 923 ,RMSE = 
27. 9283 ,P =0. 0000。 各 项 参数 值 较 模型 (2) (3) 更 好 ,所 以 ,建立 模型 (4) 

y = 0. 845361x, — 0.0120252x, + 0. 0000045x, ~ 0. 00057xs 
即 : y = 0. 845361x, - 0.0120252xs + 0. 0000045x,x, - 0. 00057xix3 (4) 
从 各 统计 量 的 数值 看 ,模型 (2) 比 模型 (1) 差 一 点 ,但 模型 (2) 比 模型 (1 ) 更 为 简单 ,模型 (3 ) 较 模型 
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(1) 2) 虽然 是 非 线性 的 回归 模型 ,但 实际 上 它 的 拟 合 效 果 并 不 好 。 模 型 (4) 是 在 加 入 交互 项 后 建立 
的 ,其 系数 中 有 的 较 小 ,几乎 为 零 , 代 人 数据 进行 预测 ,其 效果 不 是 很 好 。 综 合 上 述 情形 ,建立 的 线性 回 
归 模 型 (2) 更 能 真实 的 反映 各 自 变量 与 因 变 量 之 间 的 关系 。 


4 结 语 


本 文 从 两 方面 人 手 ,介绍 了 基于 MATLAB 的 多 元 非 线 性 回归 模型 参数 的 求解 方法 ,同时 也 尝试 了 
通过 实验 建立 多 元 非 线性 回归 模型 的 方法 ,虽然 对 于 3.2 财政 收入 预测 模型 并 没有 表现 的 像 多 元 线性 
回归 模型 那样 优秀 ,但 是 可 以 提供 一 种 思想 与 方法 ,在 以 后 的 学 习 与 实践 中 来 完善 这 种 思想 与 方法 。 
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The Multivariate Nonlinear Regression Model Based on MATLAB 


DONG Da - xiao 
(Lin Cang Teachers’ College) 


Abstract: MATLAB is the origin of the matrix operations in a highly integrated computer language. It provided a powerful scien- 
tific computing, flexible program design process, high - quality graphics and visual interface design and convenient procedures 
and language and other interface features. jn this paper, we take full advantage of MATLAB statistic tool box to achieve the multi- 
variate nonlinear regression model. Estimate of the unknown parameters. examples are used to validate the method is effective. 
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Entire Solutions of Functional Equation /5(2) +g(z) +h’(z) =1 


SU Min, Ll Yu-hua . 
( Department of Mathematics, Yunnan Normal University , Yunnan Kunming 650092 ) 


人 Abstract : This paper studies the entrie solutions of the functional equation f 6(z) + g' (2z) +h'(z) =1. Tis 
shown that there is no nonconstant entrie functions with order smaller than one satisfying the equation f(z),g 
(z),h(z) with order smaller than one satisfying the equation f°(z) +g’(z) +h:(z) =1. Meanwhile a new 
brief proof is given for the known results that if n>7 there is no nonconstant entrie functions solutions of the 
functional equation f "(z) +g"(z) +h"(z) =1. 

Key words: Functional equation; Entire function ;Value distribution theory 
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1. 学 位 论文 康 军 基于 径 癌 基 芳 数 神经 网 络 的 应 用 研究 2009 

人 工 神 经 网 络 (ArtificialNeuralNetworks， 简 称 为 ANN) 是 理论 化 的 人 脑 神经 网 络 的 模型 ， 是 模仿 大 脑 神经 网 络 结构 和 功能 而 建立 起 来 的 一 种 信息 
处 理 系 统 。 人 工 神经 网 络 的 研究 始 于 20 世 纪 40 年 代 ， 近 二 十 年 发 展 极为 迅速 ， 它 的 应 用 范围 涉及 到 数学 、 工 程 、 计 算 机 、 物 理 、 生 物 、 经 济 、 管 理 等 科 
学 领域 。 本 论文 主要 基于 ANN 的 非 线 性 逼近 性 质 ， 借 助 MATLAB 的 神经 网 络 工具 箱 ， 我 们 得 出 了 如 下 的 结果 ; 

首先 ， 我 们 实现 了 BP 网 络 〈 下 文中 简写 为 BPNN) 和 径 向 基 函 数 (RaclialBasisFunction) 神经 网 络 〈 下 文中 简写 为 RBFNN) 的 建立 、 学 习 和 训练 。 通 
过 仿真 实例 ， 体 现 了 ANN 的 非 线性 通 近 能 力 ， 并 对 BPNN 和 RBFNN 进 行 比较 研究 ， 结 果 显 示 了 RBFNN 具 有 结构 简单 、 训 练 速 度 快 ， 并 且 有 比较 强 的 抗 和 干扰 能 
， 能 达到 很 好 的 逼近 效果 。 体 现 了 RBFNN 的 函数 逼近 能 力 在 多 方面 都 优 于 BPNN。 其 次 ， 我 们 把 RBFNN 作 为 一 种 新 的 回归 方法 应 用 到 多 元 线性 回归 模型 和 多 
元 非 线性 回归 模型 中 。 通 过 实例 分 析 ， 结 果 表 明 RBFNN 方 法 用 于 回归 分 析 的 拟 合 效果 和 预测 效果 很 好 而 且 简 单方 便 。 最 后 ， 我 们 用 RBFNN 进 行 了 股票 预测 
， 通 过 以 中 国 上 证 指数 的 540 个 交易 日 的 数据 为 实验 对 象 进行 预测 ， 得 到 了 较 好 的 预测 效果 。 

本 论文 的 结构 安排 如 下 : 

第 1 章 绪 论 。 主 要 是 论文 的 研究 背景 以 及 论文 的 主要 工作 。 

第 2 章 AAN 基 本 理论 。 主 要 是 ANN 概 述 、BPNN 和 RBFNN 简 介 以 及 ANN 用 于 函数 逼近 的 一 些 定理 。 为 ANN 的 函数 逼近 能 力 的 研究 提供 了 理论 基础 。 

第 3 章 BPNN 和 RBFNN 函 数 和 逼近 仿真 。 主 要 是 利用 MATLAB 神 经 网 络 工具 箱 ， 把 BPNN 和 RBFNN 应 用 于 函数 逼近 的 设计 及 仿真 。 

第 4 章 RBFNN 模 型 与 多 元 回归 模型 。 主 要 是 RBFNN 方 法 解决 多 元 线性 回归 和 多 元 非 线 性 回归 问题 ， 并 对 传统 回归 方法 和 RBFNN 方 法 在 实际 应 用 中 进行 了 
比较 研究 。 

第 5 章 RBFNN 模 型 对 股票 市 场 价格 进行 预测 。ANN 可 以 看 成 是 一 个 具有 自学 习 功 能 的 “黑箱 ”， 在 处 理 无 法 精确 建 模 的 问题 上 具有 很 大 优势 。 由 于 股票 
的 价格 是 非 线性 的 时 间 序 列 ， 本 文 提 出 了 基于 RBFNN 的 股票 价格 预测 模型 ， 并 应 用 该 模型 对 于 中 国 上 证 指数 进行 了 预测 。 

结语 。 包 括 了 论文 的 主要 结论 与 不 足 ， 并 为 未 来 的 研究 提出 建议 。 其 中 第 3、4、5 章 是 研究 的 核心 部 分 。 


2. 期 刊 论文 汪 风 . 黄 彩 虹 . Wang Lan. Huang Caihong 基于 MATLAB 色 差 预 测 多 元 回归 模型 的 研究 -计算 机 与 应 用 化 学 
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由 于 间 和 区 染色 无 法 实现 色差 在 线 测量 , 直接 影响 染色 的 质量 . 在 染色 过 程 中 , 色差 受 多 种 工艺 参数 的 影响 , 其 中 主要 的 有 染料 的 初始 浓度 、 染 色 时 间 、 
染色 温度 以 及 助 剂 的 浓度 . 为 了 预测 色差 值 , 本 文 在 实验 数据 分 析 的 基础 上 , 利用 MATLAB 统 计 工 具 箱 的 “nlinfit” 函数 , 建立 色差 的 多 无 非 线性 回归 模型 , 并 
通过 多 元 检验 确定 该 模型 的 正确 性 . 实验 结果 表明 , 多 元 非 线性 模型 计算 得 到 的 色差 值 与 实际 值 很 接近 , 对 生产 实践 具有 一 定 的 指导 意义 . 


3. 学 位 论文 李 宗 波 船舶 转向 过 程 中 速度 下 降 问 题 的 研究 2008 

本 文 提出 的 船舶 操纵 性 指数 K 、T 估算 公式 较 原 有 的 线性 估算 公式 大 幅 提 高 了 K 、T 的 预报 精度 ， 在 保证 回归 方程 显著 性 的 基础 上 兼顾 了 数据 拟 合 
的 精度 和 回归 多 项 式 的 简洁 性 ， 符 合 精简 实用 的 原则 ， 对 操 船 者 在 航海 实践 中 评估 船舶 的 操纵 性 能 有 一 定 的 参考 价值 。 

随 厦大 型 船舶 的 数量 不 断 增 长 ， 船 舶 回转 过 程 的 速 降 问题 已 引起 了 人 们 的 重视 。 原 则 上 讲 ， 回 转 中 的 速 降 需 用 非 线 性 解析 方法 来 确定 ，Abkowitz 非 
线性 数学 模型 以 及 分 离 型 数学 模型 MMG 模 型 ) 虽然 能 预报 包括 速 降 在 内 的 各 种 回转 参数 的 变化 情况 ， 但 都 存在 计算 量 大 、 较 为 复杂 的 问题 。 因 此 ， 实 践 中 
通 间 采 用 近似 方法 来 估算 之 。 戴 维 进 、 费 加 耶 夫 斯 基 、SchoenherTr 等 国外 学 者 绘制 了 以 相对 回转 直径 为 碍 表 引 数 的 速 降 系 数 曲线 图 或 者 给 出 了 速 降 估算 
公式 ， 但 这 些 图 表 无 法 对 船舶 在 不 同 载 况 (满载 、 半 载 、 压 载 ) 下 表现 出 的 回转 速 降 差 异 给 出 合理 的 解释 。 此 外 ， 这 些 图 表 和 公式 只 能 用 于 估算 定常 回转 
阶段 速度 V 与 初始 速度 V0 之 比 ， 而 对 于 船舶 进入 定常 回转 阶段 之 前 的 过 渡 阶 段 的 速度 下 降 情 况 无 法 进行 估算 ， 而 过 渡 阶 段 恰 好 是 船舶 航行 过 程 中 进行 转 癌 
避让 经 常 涉及 的 范围 ， 因 而 以 上 学 者 的 研究 成 果 存 在 一 定 的 局 限 性 。 

本 文 另 尽 蹊 径 ， 从 由 响应 型 船舶 运动 数学 模型 推导 出 来 的 船舶 操纵 性 指数 K、T 入 手 ， 用 多 项 式 回 归 的 办 法 对 47 笨 船舶 的 K、T 指 数 统计 资料 进行 分 析 
， 得 出 了 具有 一 定 精 度 的 K 、T 指数 估算 公式 。 然 后 ， 以 此 为 基础 对 多 艘 船 的 回转 试验 速 降 资 料 进行 回归 分 析 ， 把 船舶 的 速 降 系 数 表征 为 船舶 旋回 性 指 
数 K 、 船 舶 施加 的 舵 角 以 及 航向 改变 量 人 的 函数 关系 ， 初 步 赤 服 了 国外 关于 回转 速 降 研究 成 果 的 局 限 性 ， 实 现 了 对 转 癌 全 过 程 的 速 降 情 况 的 估算 。 

本 文 以 超大 型 油轮 OPALIA 号 为 例 进 行 了 回转 速 降 的 预报 ， 用 MATLAB 绘 制 了 速 降 曲线 ， 并 与 国外 学 者 的 估算 公式 和 曲线 图 进行 了 对 比 ， 取 得 了 理想 的 
效果 ， 预 报 误差 最 大 为 14. 6% ， 最 小 为 0.3% 。 这 表明 ， 用 多 元 非 线性 回归 的 方法 实现 对 船舶 整个 转向 过 程 中 速 降 情况 的 预报 具有 一 定 的 可 行 性 。 


4. 期 刊 论文 刘 俏 . 高 伟 星 . 权 春 善 . 范 圣 第 .Liu Qiao. Gao Weixing. Quan Chunshan. Fan Shengdi 基于 MATLAB 发 酵 培 
养 基 的 优化 -计算 机 与 应 用 化 学 2006, 23 (7) 

发 酵 培养 基 配 比 的 优化 是 含有 多 变量 非 线 性 目标 函数 的 优化 . 以 青霉素 发 酵 为 例 , 通过 考察 青霉素 效 价 产 率 的 大 小 ,用 MATLAB 统 计 工 具 箱 的 
“Tegress “函数 对 培养 基 配 比 进 行 多 元 非 线 性 回归 , 以 含有 约束 的 回归 方程 作为 目标 函数 , 调用 MATLAB 优 化 工具 箱 的 "fmincon 函数 , 对 其 进行 优化 ,确定 了 
适宜 的 培养 基 配 比 . 将 优化 条 件 得 到 的 预测 值 与 实测 值 比较 , 当 培养 基 配 比 为 乳糖 8. 04、 花 生 饼 粉 4 0% 及 荣 乙 酰胺 0. 213% 时 , 目标 产物 的 效 价 产 率 为 
1226. 0, 与 优化 结果 1210. 9 比较 接近 . 


5. 学 位 论文 刘涛 水 库 水 质 参 数 了 预测 和 富 营 养 化 预警 2005 

当今 全 球 水 资源 短缺 已 经 成 为 不 争 的 事实 。 水 质 富 营养 化 问题 作为 全 球 性 的 水 污染 问题 之 一 ， 正 日 益 引起 全 社会 的 关注 。 尤 其 是 作为 城市 饮用 水 水 
源 的 水 库 富 营养 化 问题 更 是 严重 威胁 着 城市 供水 的 安全 。 本 文 以 于 桥 水 库 为 研究 对 象 ， 以 保障 城市 安全 供水 为 目标 ， 系 统 地 研究 了 水 库 富 营养 化 问题 。 

首先 ， 文 章 对 于 桥 水 库 的 水 质 现状 进行 了 调查 研究 、 收 集资 料 等 工作 。 并 在 此 基础 之 上 ， 对 所 收集 的 资料 进行 了 筛选 和 分 析 ， 绘 制 出 于 桥 水 库 水 质 
参数 的 时 空 分 布 规律 曲线 。 

其 次 ， 藻 类 恶性 增殖 作为 富 营养 化 的 核心 体现 ， 同 时 也 直接 威胁 着 城市 供水 安全 ， 所 以 本 文 以 藻类 为 切入 点 ， 选 择 反映 藻类 数量 的 叶绿素 浓度 作为 
水 质 富 营养 化 指标 和 预测 预警 对 象 进行 研究 ， 分 别 运 用 多 元 非 线 性 回归 方法 、 人 工 神经 网 络 方法 建立 了 预测 模型 分别 采 用 SAS 和 MATLAB 编 写 了 程序 进行 


















































模型 的 标定 和 验证 求解 。 进 而 采用 CE-QUALW2 理 论 模型 的 方案 加 以 辅助 验证 。 并 分 别 比较 了 三 种 方法 的 优 劣 。 

最 后 ， 针 对 上 述 研究 结果 ， 采 用 当前 流行 的 水 质 预 警 理 论 应 用 于 于 桥 水 库 的 实际 情况 ， 研 究 制定 出 了 相应 的 分 级 预警 方案 。 并 结合 2003 年 全 年 的 预 
测 数 据 给 出 预警 级 别 ， 效 果 民 好 。 

至 此 ， 本 文 结合 于 桥 水 库 的 实际 情况 系统 地 研究 了 水 库 语 营养 化 问题 ， 提 出 了 一 套 完整 的 解决 方案 ， 实 际 分 析 的 结果 显示 该 方案 切实 可 行 ， 可 以 应 
用 于 实际 。 





























6. 期 刊 论文 李 晓 燕 . 张 曼 . 余 灯 广 .Li Xiaoyan. Zhang Shu. Yu Dengguang 三 维 打印 成 形 粉 末 配 方 的 优化 设计 -机 械 
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利用 均匀 设计 方法 对 三 维 打印 成 形 的 粉末 配方 进行 了 试验 研究 , 对 试验 数据 进行 回归 分 析 , 建立 了 粉末 成 分 配方 对 制 件 性 能 的 多 元 非 线 性 回归 方程 . 设 
计 了 粉末 配方 的 多 目标 优化 模型 , 利用 MATLAB 软 件 求解 出 最 佳 配方 并 通过 试验 加 以 验证 . 采用 均匀 设计 试验 进行 回归 分 析 和 优化 是 安排 多 水 平 配方 设计 的 
良好 方法 . 


7. 学 位 论文 李 壤 水 库 水 质 再 营养 化 预测 预警 研究 2002 

水 质 富 营养 化 是 全 球 性 的 水 污染 问题 , 正 日 益 引 起 全 社会 的 关注 , 尤其 是 作为 城市 饮用 水 水 源 的 水 库 富 营养 化 问题 更 是 严重 威胁 着 城市 供水 的 安全 . 该 
文 以 于 桥 水 库 为 研究 对 象 , 以 保障 城市 安全 供水 为 目标 , 系统 地 研究 了 水 库 富 营 养 化 问题 . 首先 针对 水 库 水 质 管理 不 完善 的 问题 , 引入 了 先进 的 GIS 技 术 , 研 
制 开 发 了 相关 的 水 库 水 质 管理 软件 , 并 结合 于 桥 水 库 的 具体 地 形 图 和 水 质 信息 建立 了 于 桥 水 库 水 质 管理 信息 系统 , 实现 了 水 库 水 质 的 信息 化 管理 , 并 为 富 营 
养 化 的 预测 、 预 警 和 防治 规划 提供 了 信息 基础 . 其 次 , 党 类 恶性 增殖 作为 富 人 营养 化 的 核心 体现 , 同时 也 直接 威胁 着 城市 供水 安全 , 所 以 该 文 以 藻类 为 切入 点 
, 选择 反映 举 类 数量 的 叶绿素 浓度 作为 水 质 寅 营养 化 指标 和 预测 预警 对 象 进行 研究 , 先后 运用 多 元 非 线性 回归 方法 和 人 工 神经 网 络 方法 建立 了 预测 模型 , 分 
别 采 用 SAS 和 MATLAB 编 写 了 程序 进行 模型 的 标定 和 验证 求解 , 并 根据 预测 的 结果 结合 水 库 的 实际 情况 制定 了 相应 的 分 级 预警 方案 , 在 预警 的 基础 上 研究 相应 
的 应 急 防治 措施 , 这 些 措 施 包括 水 库 内 的 党 类 前 减 措施 和 水 三 采取 的 工艺 和 药剂 方面 的 应 急 准 备 . 最 后 , 为 了 真正 实现 安全 供水 , 对 于 桥 水 库 富 营 养 化 成 因 
进行 了 分 析 , 并 对 现 阶 段 可 以 实现 的 防治 措施 进行 了 系统 的 规划 , 使 用 Lingo 软 件 对 防治 规划 进行 了 优化 分 析 , 在 保证 水 质 、 满 足 环境 容量 约束 的 条 件 下 得 
到 了 最 优 的 富 营 养 化 防治 规划 . 
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